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The most general covariant action describing gravity coupled to a scalar field with only second 

order equations of motion, Horndeski's theory (also known as "Generalized Galileons"), provides an 

all-encompassing model in which single scalar dark energy models may be constrained. However, the 

fT^ generality of the model makes it cumbersome to manipulate. In this paper, we demonstrate that 

T-H when considering linear perturbations about a Friedmann-Robertson- Walker background, the theory 

^^ is completely specified by only six functions of time, two of which are constrained by the background 

O^ evolution. We utilise the ideas of the Effective Field Theory of Inflation/Dark Energy to explicitly 

K^,^ construct these six functions of time in terms of the free functions appearing in Horndeski's theory. 

rrt These results are used to investigate the behavior of the theory in the quasistatic approximation. 

k^ We find that only four functions of time are required to completely specify the linear behavior of 

^—{ the theory in this limit, which can further be reduced if the background evolution is fixed. This 

/-jv presents a significantly reduced parameter space from the original presentation of Horndeski's theory, 

I giving hope to the possibility of constraining the parameter space. This work provides a cross-check 

for previous work on linear perturbations in this theory, and also generalizes it to include spatial 

I I curvature. 
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a I. INTRODUCTION 

\^ Since the discovery of the accelerated expansion of the universe (TJ [5] , a plethora of models to explain this 

J^ phenomena have been proposed (see Ref. |3J for a comprehensive review). In order to obtain a homogeneous 

I and isotropic expansion, the mechanism behind the accelerated expansion is typically attributed to a (effective) 

r- — . scalar field. It is, however, generally straightforward to select parameters for a model to select the desired 

\^ expansion history of the universe [4j, and so it is important to investigate the perturbative behavior of a 

■ ' model in order to establish observational limits on the parameter space. Unfortunately, this analysis is 

f— ^ typically much harder to perform than the background analysis. 

CO Recently, there has been much interest in general models that attempt to contain a number of individual 

'""' models [SHU] • One of the motivations for these approaches is to be able to construct model- independent 

^ constraints, which will hopefully save time and effort analyzing individual models. This paper addresses the 

intersection of two such methods. 

An old result by Horndeski [T^], recently rediscovered by Deffayet et al. [HKHj, states that the most 
C^ general possible covariant scalar field theory coupled to gravity without higher order derivatives is given by 

the action 



'^ = E^* = E / d'-T^^^^- (1) 
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- Gs.jf [{a^f - Hacp) (v^v^ifi) (v^v^^) + 2(v^v„^) (v"v^</.) (v^v^./.)] . (2d) 



where 

C2 = K{ip,X), (2a) 

C3^-G3{ip,X)nip, (2b) 

£4 = Giiip, X)R + G4,x [(0(^)2 - (V^V^ip) (V'V^)], (2c) 

/:5 = G5(^,X)G^,(V''VV) 
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Here, X = —g^"d^Lpd^ip/T\ and K^ G3, G4 and G5 are arbitrary free functions^ This action wiU henceforth 
be referred to as the "Horndeski action" . The generahty of this action suggests that almost all single scalar 
field dark energy models are contained within it (Lorentz violating models such as Hofava-Lifshitz gravity 
[Is] are the biggest category of exceptions). On the other hand, this description is almost too general, in 
that the impressive amount of freedom contained in the action prohibits the calculation of observational 
constraints for the model. 

A particularly beautiful approach to perturbations in cosmological spacetimes comes from the Effective 
Field Theory (EFT) of Inflation [16j . In this approach, it is assumed that the background evolution of the 
universe is known, allowing the perturbative behavior to be investigated in detail. This approach has been 
generalized to dark energy models [S| 'T] TT] , where it has been demonstrated to provide a very general method 
for calculating perturbative behavior. By discarding information about the background, the EFT approach 
condenses the model-dependent information into a handful of functions of time. It is hoped that by placing 
model-independent constraints on these functions of time, we will be able to constrain very general theories 
of dark energy, such as Horndeski's theory. 

The goal of this paper is to perform a linear perturbation analysis of Horndeski's theory about a Friedmann- 
Robertson- Walker (FRW) background by casting it into the formalism of the EFT approach. De Felice et al. 
[18] have already performed such a feat by calculating the linearized equations of motion and investigating the 
results. Here, we simplify matters somewhat by demonstrating that Horndeski's theory can be condensed from 
its original four functions of two variables down to just six functions of time, for the purposes of describing 
linear perturbations. Furthermore, we generalize upon the results of De Felice et al. by performing the 
analysis with arbitrary spatial curvature. Once the results have been obtained, we use the machinery already 
developed for the EFT approach to calculate the modifications from general relativity in the quasistatic limit. 

This paper is arranged as follows. We begin by providing a very brief overview of the EFT of Inflation 
formalism in Section [Hj Next, we use counting arguments to decide which operators will be necessary to 
describe Horndeski's theory within the EFT framework in Section |III[ We then demonstrate how the two 
approaches can be matched together, and explicitly detail the equivalence in Section [TV] Finally, we investigate 
the perturbative behavior of the theory in the quasistatic approximation in Section|V{ paying close attention 
to the number of free functions required for a complete description. We conclude by discussing the application 
of these results to inflationary models and suggest future directions. 



II. OVERVIEW OF EFFECTIVE FIELD THEORY APPROACH 

The EFT of Inflation formalism [TS] provides a description of perturbations in single scalar field theories 
about an FRW background. It works by using the unitary gauge, in which the foliation of spacetime is 
chosen such that the scalar field has no perturbations (the scalar degree of freedom has been "eaten" by 
the metric) . Because of the symmetries of FRW, the possible operators in an action in this gauge are very 
straightforward, and can be arranged into a perturbative expansion, with each operator having a function of 
time for a coefficient. In this section, we provide a very brief overview of the formalism, directing the reader 
to the literature for a thorough introduction [5J [71 [TS] . 



^ Note that there exist different definitions for X in the literature, variously including/excluding the factor of 1/2 and the minus 

sign. 
^ We omit a matter action, which we take to be minimally coupled to the metric in the Jordan frame. This assumes that the 

weak equivalence principle holds. 



Begin by assuming that the background metric takes the form 



(3) 



where gij is a maximally symmetric time-independent spatial metric with scalar curvature R^^' — 6kQ. The 
operators that can be included in the unitary gauge action must be constructed of the following objects. 
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These are, in order, functions of time, a perturbation to the 17™ component of the inverse metric, the spatial 
metric, a perturbation to the extrinsic curvature tensor to surfaces of constant tima^ a perturbation to 
the Riemann tensor of the spatial metric, covariant derivatives associated with the spatial metric, time 
derivative^ and the three-dimensional antisymmetric tensor. Evaluated on the background, the spatial 
metric is given by hij = a^gij. Because there are no e symbols appearing in Horndeski's theory, they will 
not be required in the perturbation description, and we ignore them from now on. The perturbations are 
constructed by subtracting the background values written in terms of just functions of time, the spatial 
metric and the normal vector. 

It was previously shown [5] that the leading order operators in such a construction are given by the following 
action. 
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The functions ri(i), A(i) and c{t) are related to the background evolution of the comology, and appear in the 
Friedmann equations as [HJ [7] 
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where pm and Pm are the energy density and pressure arising from matter and radiation fields, H is the 
Hubble parameter d/a, and overdots refer to time derivatives. 

Through the use of a technique known as the "Stiickelberg trick" , the scalar field perturbations can be 
reintroduced to the action (Is]). The perturbations are typically called tt, which relate to the scalar field as 

(p = ^°{t + Tr) = ip°{t) + S^ (8) 

where (^" is the background scalar field. Both the equations of motion for the tt field and the stress-energy 
tensor were calculated in [6]. 



III. CHOICE OF OPERATORS 

In order to construct linear perturbations in Horndeski's theory within the EFT formalism, we first need to 
identify which operators will be required. Using the notation in the action (l5|, the operators with coefficients 



Different conventions exist for the sign of the extrinsic curvature tensor. We use the definition K/^i, = {5^ + n^n^)VA»^iy, and 



define the normal to surfaces of constant time as riu 



^ d^t/^- 



nOO 



C^, where C^ is the Lie derivative along the shift vector 



^ Technically, the object that is invariant under the symmetries is dt 
in an ADM decomposition. However, because the shift vector is vanishing on the FRW background, the terms involving the 
shift are higher order in perturbations and can be ignored. 



0(i), A(i) and c(i) will definitely be required, as these operators specify the background FRW evolution. The 
remaining operators, which only affect the behavior of the perturbations, should have the property that they 
contain only two derivativefP] Of course, we may take linear combinations of the remaining operators to 
cancel out any higher order derivatives. 

For linear perturbation theory, we are interested in operators that are second order in perturbed quantities. 
Furthermore, all of these operators must be constructed out of the objects listed in Eq. Q. Consider the 
metric in Newtonian gauge, which contains no derivatives acting on the scalar perturbations. The object Sg'^^ 
has no derivatives, nor does the spatial metric, while the Riemann tensor contains two spatial derivatives, and 
the extrinsic curvature tensor (being the Lie derivative of the spatial metric along the normal to hypersurfaces 
of constant time) contains one time derivative. Therefore, we can join these objects together to make a 
complete list of all possible operators at second order in perturbations that contain two derivatives or lesfrl 

iSg^r, 5g"'SK\, {5K\)\ 6K)K^,, Sg°°6R^'\ D,Sg°'D%'', {d.Sg^'f, d,Sg''SK\ (9) 

Here, indices have been raised and lowered using the spatial metric, and we have taken advantage of integration 
by parts. 

Unfortunately, the "two-derivative" property is a little muddled, because once the Stiickelberg trick has 
been applied, the number of derivatives acting on the Stiickelberg field can be different from the rest of 
an operator. Here is how the Stiickelberg field is introduced to each object, truncated at first order in 
perturbations. 

<5(7"° ^ <5g°° - 27T (10) 
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In Table IT] we count the number and type of derivatives acting on the metric and scalar field perturbations 
for each of the operators in ^. 
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Table I. The number of space and time derivatives in diflterent operators, sorted by whether the derivatives are acting 
on metric perturbations only, scalar perturbations only, or a mix of the two. Only the highest order derivatives in 
each category are listed. 

- No derivatives are present 

ns There are n spatial derivatives 

nt There are n time derivatives 

-I- Multiple combinations of derivatives are present at the highest order 



^ At least, at second order in perturbations. At higher orders, more terms can become involved in order to provide the required 

cancellations. 
® Gleyzes et al. |19) point out that this condition is actually too restrictive; it is possible for a combination of operators containing 

higher derivatives to have their higher derivative terms cancel. They show that the combination W '^■' 5K^- — i?,('''<5A'V/2 can 
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is a total derivative, even 



be reduced to operators in this list however, and that the combination 2(R"')^ — 3/4iJ 
when multiplied by a function of time (note that in both these cases, the curvature scalar and tensor are represented without 
the background component subtracted off). We are unaware of further examples, and it is likely that in any case, such may be 
rewritten (as these may be) in terms of the operators presented here. 



We see from the table that the operator {Sg°°)'^ is fine, as only second order derivatives are present. The 
operator Sg'^^5K\ might be problematic because of the cross term —2'kD'^n, but the time derivative can be 
removed using integration by parts, and so this operator is also fine. 

The operator {dtSg'^^)^ is alone in that it contains four time derivatives acting on scalar perturbations, 
so these derivatives cannot be cancelled in combination with another operator, and this operator must be 
dropped. Similarly, the operator dt5g^^5K contains terms with three time derivatives, which cannot be 
cancelled with anything else, and so this operator must also be discounted. The operator DtSg^^D'^dg^'^ 
contains terms with two time derivatives and two spatial derivatives which cannot be cancelled in combination 
with any remaining terms, and so this term is also thrown out. 

The operators [SKY and SK'^^SK-'- both contain terms with four spatial derivatives acting on the scalar 
perturbations. A simple analysis reveals that so long as the coefficients of these operators are equal but 
negative of each other, these terms cancel. Both operators also contain terms involving one time and 
two spatial derivatives. Looking at the table, it is feasible for the operator Sg^^SR'^^^ to provide an exact 
cancellation for this term. Indeed, these terms do cancel when the coefficient of Sg'^^SR'^'^^ is given by 

2Af2 = Mi = -Af|. (13) 

To see this cancellation at the quadratic order requires expanding the three operators in metric perturbations 
and using integration by parts, which we demonstrate in the following section. The operator Sg'^^SR^^' also 
has a cross term ttD^tt which can be dealt with through integration by parts, just as for Sg^^'^dK^^^. 

Thus, we expect that linear perturbations in Horndeski's general theory about an FRW background will be 
described in the EFT of Dark Energy context by the following action. 

s = l d^xy^l "fnm + m - <tw + ^ (sg'r 

M!{t),nOs:r.^ Af|(t) 



SgOOsK\ - ^ (^SK\^ - 5K^^5K\ + 2Sg"^5R^'^) [ + S^atteAg,.] (14) 



We claimed this result without proof in [6] . All that now remains is to identify the correspondence between 
the six functions of time appearing in this action and the free functions appearing in the Horndeski action. 



IV. MATCHING HORNDESKI TO THE EFT CONSTRUCTION 

The Horndeski action is sufficiently complicated that performing explicit matchings between the various 
functions by hand is a daunting task. Instead, we pursue a different approach, turning to computer algebra 
software. 

As we are considering linear perturbations about a background in Horndeski's theory, we need to consider 
quantities that are quadratic in perturbations in the action. It is thus sufficient to expand the action for 
Horndeski's theory to quadratic order about a background, and compare this to the EFT action expanded to 
the same order. By matching coefficients in the action, we can thus express the functions of time in the EFT 
formalism in terms of background quantities in Horndeski's theory. Note that the leading order piece of this 
expansion (zeroth order in perturbations) is a meaningless boundary term which can be discarded. 

We work in unitary gauge in which there are no scalar field perturbations, and so the scalar field is given 
hy (fi = ip^{t) (the superscript zero will be dropped from now on). This gauge requires us to consider a metric 
with three scalar perturbations, which we write as 

ds^ = -(1 + 2(j))dt^ + d.Bdx'dt + a'^g^J{l - 2V')da;Mx^ (15) 

The three scalar perturbations are 0, ip and B. 

When we expand the action in terms of these perturbations, there will be a number of different terms that 
are quadratic in the perturbations. However, some of these can be manipulated using integration by parts. 
We will preferentially integrate by parts to move spatial derivatives onto B where possible, then ip. Similarly, 
we will move all time derivatives onto ^. 



A. Quadratic Expansion: EFT 



The EFT action we wish to expand is 
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where we have left the last three functions of time separate in order to demonstrate the required cancellations 



of higher derivative terms. Using the metric (15), it is straightforward to express the perturbed objects in 
this action to the appropriate order. 
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The expansion of the Ricci scalar, including the metric determinant and integrating by parts, is a straightfor- 
ward but somewhat tedious calculation. 
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Here, R° ^ 6H + 12H'^ + efco/a^ is the background Ricci scalar, and D = {l/a'^)g'W,Dj. 

Combining all of the terms in the action yields the following expression for the linear and quadratic terms. 
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Note that although there are thirteen coefHcients to match at quadratic order, only eleven of these are 
independent. As we have eight functions of time in the action, this leaves three extra equations to ensure 
consistency. 

Two operators in this action should look problematic from the perspective of derivative counting. The most 
obvious is the last, DBDB, which evidently has too many derivatives, enforcing M2 = —M^. The other is 
ipDB. Here, however, the precise issue of how many derivatives are present is clouded by the gauge. Once the 
Stiickelberg trick has been performed (which changes B^tt, ip^(f>— TTjip^ip — Hit, in Newtonian gauge) 
the term (j)^\^jJ yields a term ttDt/), which combines with ipOB. The requirement that this three-derivative 
operator cancels yields 

4Af 2 - 3M| - Aff = 0, (23) 

giving the result 2Af ^ — Afl stated previously. 



B. Quadratic Expansion: Horndeski 

We investigate each of the terms in the Horndeski action in turn to demonstrate how they can be written in 
terms of operators in the EFT formalism. Note that the contributions from different terms to the coefficients 
of each operator in the EFT add linearly, so it is sufficient to investigate them in isolation. 

L2 

We begin with the term £2 = K{ip,X). The mapping of this operator to the EFT construction has been 
demonstrated in the context of fc-essence models in a variety of papers, including [TllinillZ]- Because the 
construction is trivial, we do not use the machinery from above. 

In unitary gauge, X ~ {1 — Sg°°)Xo, with Xq = ip'^/2. Therefore, we simply expand K in powers of Sg'^°, 
using the chain rule as necessary. 

K{4>, X) = Ki^, Xo) - K^xX^Sg"" + ]^K,xxXl(5g''''f + . . . (24) 

Here, Kx indicates a derivative of K with respect to X (such derivatives are always evaluated on the 
background ((/),Xo)). We can then read off the contribution to the following coefficients, with all others 
vanishing. 
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We next look at the term £3 = —03(1,5, XyOnp. This term has been previously cast in the EFT language in 
[7], although only results under a simplifying assumption were presented. The expansion of this term (and 



the following two) into the form of ( 22 1 is rather long and uninsightful, so we present only the coefficients 



that we have extracted by matching the results with the EFT expansion 
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c(t) = -^(^G3,x(3i/^' + w)-2G3j (26) 

A(t) = ~G30 (27) 

Ml{t) = J^2 (^G3,x(3i/^ + <^) + G3,xx^'(3i?95 + (/3) - G^^x^) (28) 

Ml{t) = -G^^xV^ (29) 

Mlit) = -Miit) = 2M^{t) = n{t) = (30) 

Note that we choose to eliminate partial derivatives with respect to if in favour of time derivatives, where we 
use the chain rule 

G3 = G3,^(y9 + G3,xW- (31) 



To our knowledge, this term and the fohowing have not previously been cast in the language of the EFT 
formalism. This is the first term that requires all operators to be present. 
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C. Summary 



We have demonstrated that a combination of EFT operators, chosen such that at most second order 
derivatives exist in the action, is capable of describing linear perturbations to the Horndeski action about an 
FRW background. The matching was performed very generally, including spatial curvature and an arbitrary 
background scalar field p{t). Although ip>{t) must satisfy ip > for the unitary gauge to exist, we conjecture 
that the resulting EFT description (once transformed out of unitary gauge) will match Horndeski's theory 
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Table II. The EFT operators needed to match each term in the Horndeski action. 

even when this condition does not hold. The EFT operators required to match each term in the Horndeski 
action are summarised in Table llll 

Although we have only performed the matching using scalar perturbations, given that there are no other 
operators with at most two derivatives, the matching must also hold for tensor and vector perturbations. 
Thanks to the decomposition theorem, the scalar, vector and tensor modes are decoupled at linear order. 

V. OBSERVABLES 



Given the action (14), we may proceed to transform out of unitary gauge, and calculate the equations 



of motion. The procedure for the transformation is detailed in Refs. [SI [71 ITC]. while the resulting scalar 
equations of motion are presented in [6|. Here, we use these results to derive the modified Poisson and 
anisotropic shear stress equations, as well as the scalar field equations of motion. 

We work in Newtonian gauge, using the metric 



ds^ = -(1 + 2(/-)dt2 + a^gij{l - 2i!)dx^dx^. 



(44) 



We follow the conventions of Kodama and Sasaki [20] for the spatial mode functions in curved space. For 
each fc, define Yj:{x^) to be a solution of the equation 

g'W.D.Yj: = -k'Yj: , (45) 

where the spatial covariant derivatives are now those associated with the spatial metric cjij (these derivatives 
are the same as for the spatial metric hij previously, when evaluated on the background). We will typically 
suppress the k dependence of Y . Both Y and Y* will be solutions of this equation, as will any linear 
combination. We fix this freedom by choosing 



lim y(xM = e^ 
feo-s-o 



(46) 



so that the modes become the usual Fourier modes when the background spatial metric is flat. The 
normalization of Y is chosen such that 



jd^x^gY^Yl 



{2Trf S^{k-k'). 



Taking derivatives of Y , vector and tensor mode functions are defined as 



Y, = -k-^D,Y, 
Y, = k-^D,,D,Y 



1 



'-,9tjY- 



(47) 

(48) 
(49) 



We raise and lower indices on Yi and Yij with g.^ only. This ensures that Y, Yi, and Yij are independent of 
time, no matter the position of their indices. 

We use these mode functions to decompose the matter stress-energy tensor in momentum space as 

T%^~ipra+SpY) (50) 

T° = (p,„ + P„,)vY, (51) 



n = (Pm + SPY) 5) + P^UY) 



(52) 
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We define 



O TT 



as the rest-frame density perturbation of matter. 



Using the results of [B], the modified Poisson equation arising from the action (14 1 is given by 



rriQ^ 



2k^ - 6ko 



ip = — Pr?j A — 2c(7r — 



mnCl 



1/2 

-3HTr+—Tr-3('iP + H(l> 



-4M|(^- 



fc2 

Aff ( 3xp + iHcj) + 3i?7r - — tt 



2M 



2 3fco-fc^ 

2 2~ 



-V. 



(53) 



(54) 



Similarly, the modified anisotropic shear stress equation is given by 



k'^ ■ k"^ - k"^ I r 



Mo 



TT + 4> 



(55) 



Finally, the tt equation of motion is given by 



= CTT — C(/) + C — TT + C (tt - 



a' 



ScHtt — Scijj — 6cH(j) 



4 / 2 



-6(2i/2 + H)(P - 3H^ + 



UH^ -i^-2 



{k^ - 3fco) 



v- 



2Af2^ 






M| 



,A/2 



^-'^+i^a|+^^M^ 



3^ (3 fv- + H(f>) + 3Htt - —TT 

1,2 )l2 

+ 6ii0 + 3iI0-i7— TT-— , 



9iJ fvJ + H(j)] + 3(H + 3HH)tt + S-ip 



6 H- 



(V' + i?'/') 



,vfc' 



-4iJ- 



GH-" 



M^ 



{k"^ ~ 3/co) 



^+^tr^+r'+'^S+^r 



(56) 



A. Quasistatic Approximation 

We now investigate the equations of motion in the limit of the quasistatic approximation, following Refs. 
[51 [TB]. The quasistatic approximation typically consists of two approximations, namely that time derivatives 
of perturbations are small compared spatial derivatives, and also the sub-horizon approximation k/aH ^ 1. 
We refer the reader to Ref. [9 for a detailed discussion on these approximations. 

Under these approximations, the Poisson equation yields 



{2mln + 2A-/|) i) - (mlil + Mfj tt 
the anisotropic shear stress equation simplifies to 



P^"^' 



(57) 



-^. 



Mj 



Mj 



H 



2M2 



0, 



(58) 
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and the tt equation of motion becomes 






= 




^ + 


ml- Ml 
2 2 


4' 


+ 


2 y Ml J \ a^ A/2 


- 


f 


HnR0-3cH+'f 


/ 
\ 


M \ 

-J- + if + 3HH - 3M 

Ml J 




2H- 



,M2 

M. 



H + 3M|iJ2 



fc2' 



(59) 



We leave in a mass term for the tt field, as this may be comparable to k^ ja} . These equations can be combined 
in a matrix equation, 



A^ 



A. 



L-i/i L/0 L'TT + '}^^-k2 



( ^ 



V- 



^/o A 

^2 Pm'-^ 






(60) 



where A-j; are the coefheients in the Poisson equation, B^ are the coefhcients of the anisotropic shear stress 
equation, and C^ are the coefheients in the tt equation of motion. Denoting this matrix as M. , the equation 
can be inverted to obtain 



^ 




(61) 



We can now compute the effective Newtonian constant and anisotropic shear stress phenomenological 
functions, defined by 



^47rGcff(fc,t)pp„,A, 



i/) = 7(fc, i)0 . 



(62) 
(63) 



These functions have been given a variety of names by different authors, and have been investigated in a 
number of papers (see, e.g., [S1[H1[52]). For Horndeski's theory, these functions are given by 



Gcff(fc,t) 



47r 



\^~%. 



-DttG^ 



B^Ct^ 



-S^A07r2"p" 



4"" A^ [B^C-^ + B^Ct^2^ — BttC^J + Aji- {BjpC^ — B^Cjj;) 






B^Ct: — Bt^C^ + B^CTr2]^ 
B^C^ — B^jCt^ — -BifcO,r2ir 



(64) 



(65) 



fe2 



As stressed by Silvestri et al. [9 , the numerator of Gcfr and the denominator of 7 are the same. Both 
numerator and denominator are polynomials in a/fc, and as expected for a scalar field theory, only even 
powers of a/k appear. Furthermore, as expected from Horndeski's theory, only a constant and the second 
power of a/k are present. Finally, we see that Goff and 7 can be written in terms of six functions of time, in 
the form 



„ _ 1 /l +./2f2 

^-^ h + fiw 



7 = 



f5 + hj: 

h + /2I 



fc2 



(66) 



fc2 



This can be reduced to only five background functions of time by defining gi = fi/ fj for a given j. 

Comparing the function counting, we showed above that the construction of the Horndeski action in the 
EFT formalism requires six free functions of time. However, note that the function A(t) does not appear in 
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any of the equations of motion (l54|, (55) or (56 1. Furthermore, in the quasistatic hmit, the function M2 does 



not appear in the matrix A^, and therefore only four free functions are required to specify GcS and 7 for 
Horndeski's theory (assuming H is given). 

If the background evolution is completely specified a priori, the Friedmann equations (|6| and ([7]) can (in 
principle) be used to eliminate two of r2(i), A(t) and c{t) in favour of known background functions of time. 
Choosing these to be fl(t) and c(t), this implies that only two free functions of time are required to specify 
the behavior of Horndeski's theory in the quasistatic limit. This is a significant reduction in theory space, 
and bears further investigation. However, in practice, it is easiest to eliminate c{t) and A(t), leaving three 
free functions of time. 



VI. CONCLUSIONS 

In this paper, we sought to describe linear perturbations in Horndeski's general scalar field theory about 
an FRW background, using the machinery of the EFT of inflation. We began by motivating the terms in the 
EFT action that would be required by using a derivative-counting argument, and suggested that only six 
functions of time would be needed. We then used a perturbative expansion to explicitly match the coefficients 
of the EFT operators to the free functions in Horndeski's theory. 

Once the theory was constructed, we presented the modified Poisson equation, the modified anisotropic 
shear stress equation, and the scalar field equation of motion. These results provide an independent check of 
the work of De Felice et al. 118] . as well their generalization to include spatial curvature. 

We next investigated the behavior of the model in the quasistatic approximation, and used this limit to 
calculate the effective Newtonian constant and the ratio between the gravitational scalars tp/^k- It was found 
that in this limit, only four of the functions of time in the EFT construction appeared in the equations. This 
presents a reduction beyond the function counting that Silvestri et al. achieved from general arguments, 
although their arguments are more generally applicable. 

It was discussed that by fixing the background evolution of the cosmology a priori, further reductions in 
the number of functions of time could be achieved, implying that even fewer functions of time are required 
to specify the linear behavior of the most general scalar field theory. This motivates investigating how a 
principle component analysis could best constrain this theory over this reduced parameter space, which is a 
significant reduction from four functions of two variables. 

Although these results have been directed towards dark energy models, the formalism applies equally well 
to inflationary models (indeed; this is where the formalism was originally applied, with much success). There 
are two major differences when applying the formalism to inflationary models. Firstly, there is no matter 
present, and so the matter action can be ignored. Secondly, one can perform a conformal transformation 
of the metric to set 0(t) — > 1. The Friedmann equations then uniquely specify A(t) and c{t) in terms of H. 
Thus, linear perturbations to Horndeski's theory applied to inflation are described by just four functions of 
time: H, M2, Mi and M2- This is further reduced if H(t) is specified a priori. Unfortunately, inflationary 
models are typically more interested in extracting the three-point function, which go beyond our present 
results. Attempting a matching in the action at cubic order is an incredibly daunting task, although we 
suggest that the terms present at cubic order in the EFT expansion should be able to be determined by using 
similar arguments to the derivative counting arguments presented here. 

We are currently using the EFT approach discussed in this paper to investigate the parameter space 
of Horndeski's theory, and compare it to observational data, using CAMB [53^ to numerically evolve the 
equations of motion. We hope that data from future experiments such as the LSST, the Dark Energy Survey 
and Euclid will allow us to place stringent constraints on even the most general of theories. 
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Addendum: While this work was in its final stages of preparation, |19j was published, presenting similar 
results but using a different analysis. We have checked that all results are in agreement. 
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